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^ ■ Abstract 



A family of degenerate domain wall configurations, partially preserving supersym- 
metry, is discussed in a generalized Wess-Zumino model with two scalar superfields. 
We establish some general features inherent to the models with continuously degen- 
erate domain walls. For instance, for purely real trajectories additional "integrals of 
motion" exist. The solution for the profile of the scalar fields for any wall belonging 
to the family is found in quadratures for arbitrary ratio of the coupling constants. 
For a special value of this ratio the solution family is obtained explicitly, in terms of 
elementary functions. We also discuss the threshold amplitudes for multiparticle pro- 
duction generated by these solutions. New unexpected nullifications of the threshold 
amplitudes are found. 



1 Introduction 

Existence of several degenerate supersymmetric vacua is a generic phenomenon in supersym- 
metric theories with scalar superfields. Moreover, in many instances the vacuum manifold 
consists of several isolated points. Thus, the possibility arises of domain wall configurations 
interpolating between these vacua. 

It has been recently showni I that some of the domain wall configurations in (3+1) 
dimensional theories posses distinct supersymmetric properties: 



i) they generate a central extension of the N = 1 superalgebra. The wall tension is 
proportional to the central charge. Due to the non-renormalization theorem for the central 
charge this implies that the wall energy density is exactly calculable, it is not renormalized 
by the loop corrections. 

ii) they preserve two out of four original supercharges ("N=l/2 super symmetry" ) corre- 
sponding to the minimal supersymmetry in the (2+1) dimensional space tangential to the 
wall; 

Hi) the profile of the fields across these walls satisfies first order differential equations 
analogous to the Bogomolny-Prasad-Sommerfeld (BPS) equations^! (the so called "BPS- 
saturated walls", or "BPS walls"). These equations were called the creek equations® be- 
cause of a mechanical interpretation [] 



A non- vanishing central charge of the N = 1 superalgebra exists for any field configuration 
interpolating between two distinct vacuad. Not every such configuration is BPS saturated, 
however. Those domain walls that are BPS saturated posses peculiar features following from 
the point ii above. The BPS domain walls have interesting properties even in the simplest 
theories, e.g. in the Wess-Zumino model with one chiral superfield. Even more remarkable 
they become in the theories with two or more chiral superfields interacting with each other. 

In the theories with K scalar superfields generically there are 2 K degenerate vacuum 
states and, correspondingly, there are at least 2 K ~ 1 (2 K — 1) domain wall types@. In Ref. 
TOfl it was shown that, quite typically, some of these domain walls turn out to be continuously 



degenerate. Collective coordinates exist corresponding to a continuous deformation of the 
internal structure of the wall. Varying these coordinates we change the wall structure leaving 
the wall energy density intact. 



Previously the creek equations were considered in different contexts in Refs. 



1 



In this work we investigate this phenomenon, both in general aspect and in some simple 
examples. It will be argued that the continuous degeneracy is related to the existence of 
additional "integrals of motion" . We will explicitly find such an integral in a particular two- 
field model. Using this result it turns out possible to obtain a generic family of the domain 
wall solutions in this two-field model. For arbitrary values of the coupling constants the 
solution is given in terms of quadratures. For some specific values a closed-form solution in 
terms of elementary functions exists. We take advantage of this explicit solution to extract 
consequences for the high-order multiparticle amplitudes at threshold. 

The supersymmetric vacua are determined by the extrema of the superpotential W(<ftk), 

dW/dfa = 0, k = 1...K. (1) 

The general form of the BPS-saturation equations for a static wall in which the fields (pi 
depend only on the coordinate z is 

dW ] ■ 

dAk = — f e M , (2) 

where a is a constant (z independent) complex phase. Let us assume that two solutions of 
the equations (|l|) are found, {4>}i and {4>} 2 , where the braces denote a set of all scalar fields 
(p k . Denote the corresponding values of the superpotential by 

Wi = W{{(j)}i) , W 2 = W({<f)} 2 ) . (3) 

Without loss of generality one can assume that W± = and W 2 is real and positive (this can 
be always achieved by appropriate transformations of the superpotential). Then the phase a 
in Eq. (0) can be set equal to zero. If, additionally, the superpotential is real for real values 
of the superfields, to which we limit our investigation, then Eqs. (0) take the form 

dW 

dAu = ~~prr~ • (4) 

Now, if one interprets z as "time", the latter equations have a simple mechanical inter- 
pretation@: they describe the flow of a very viscous fluid, whose inertia can be neglected, in 
the potential relief — W(4>k) from one extremum — W\ of the relief, along a gradient line, to 
a lower extremum —W 2 . (Obviously, reflection of the z direction is possible, in which case 
the flow is described by the "potential" W rather than —W. Instead of the reversal of the 
z axis one can view this as a shift by it of the phase a in Eq. (^j. For definiteness we use 
the mechanical analogy within the conventions of Eq. (p.) From the fluid flow analogy it is 
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clear that the necessary condition for the existence of a solution of the equations (£|) is that 
the extrema W\ and W% be of different height. One can take advantage of a rich intuition 
one has in connection with the mechanical motion of this type. To make this analogy more 
graphic we will sometimes denote the derivative over z as 

d z (f) -> 

in the cases where there is no menace of confusion. Correspondingly, the quantities conserved 
along the given trajectory will be referred to as integrals of motion. One integral of motion, 
energy, is well-known of course; it is universal and has nothing to do with the specific 
trajectories under consideration. We will be interested in additional integrals of motion, 
specific to the creek equations. 

The surface energy density of the BPS wall is given by the difference of the superpotential 
at the two extrema, 

e = 2(W 2 -W 1 ) . (5) 

Although at least some of the results to be reported below seem to be generic, we will 
make no attempt at formulating them in full generality. Instead, we will focus on a specific 
instructive example: a two-field modelB E9 with the scalar superfields <3> and X and the 
superpotential 

W{&, X) = — $--A$ 3 -a$X 2 . (6) 
A 3 

Here m is a mass parameter and A and a are dimensionless coupling constants. It is assumed 

that the phases of the fields and W are adjusted in such a way that all parameters in the 

superpotential are real and positive. Only occasionally we will digress to more general 

models. Below the lowest component of $ and X will be denoted by and \i respectively. 

Clearly, the superpotential (j^) is not the most general form of the superpotential in the 

renormalizable (3+1) dimensional models, even given the freedom of redefining the fields. 

It contains, however, sufficient features for a discussion of the non-trivial properties of the 

domain walls we are interested in. Q The model with the superpotential (Q) will be referred 

to below as a minimal two-field model. In terms of the scalar components and x °f 

the respective superfields $ and X the "potential relief — W has its maximum — W\ at 

= — m/A, x = and the minimum —W2 at = m/A, x = 0. It also has two saddle 

points of equal height at = 0, x = A a. The BPS-saturated walls exist connecting 

2 An equivalent model was also considered in Ref. jllj , and a similar model of higher order in the fields 
in Ref. ff^|, in connection with different problems. 
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the maximum and the minimum, and also connecting either of the saddle points with the 
maximum or the minimum. Moreover, all these BPS-saturated configurations belong to 
one and the same family of solutions, corresponding to the flow from the maximum to the 
minimum with different starting conditions©. All trajectories from the family are real. 

This is an ideal setting for establishing the existence of additional integrals of motion. The 
one relevant to the model (||) is obtained in an explicit form. Certainly, given the additional 
constraint, one can readily reconstruct the full family of solutions, in quadratures. Further 
simplifications arise (a) for arbitrary values of A, a and the vanishing value of the integral of 
motion; (b) for A/a = 4 and arbitrary value of the integral of motion. As a matter of fact, 
the first case was treated in Ref. |I| where it was found 




Up till now, this was the only non-trivial solution for one specific configuration in the de- 
generate family that was explicitly known, apart from the trivial standard wall with \ = 0. 
In this work in Sect. 3 we construct the solutions for all configurations belonging to this 
family. For arbitrary ratio of the coupling constants p = X/a this solution is semi-explicit 
in the sense that the trajectory in the field space is found in terms of elementary functions, 
while the dependence of the fields on z, although expressed in quadratures, is not found 
analytically. An explicit dependence on z can be found in terms of elementary functions for 
two special values of p: p = 1 and p = 4. The case of p = 1 is, however, trivial since after a 
7r/4 rotation in the space of the fields (<p, x) the model reduces to two fields not interacting 
with each other. The remaining case of p = 4 is quite non-trivial and provides a whole family 
of new domain wall solutions. 

In four dimensions the requirement of renormalizability restricts the form of the super- 
potential: it must be polynomial in fields of at most third order. If we dimensionally reduce 
the theory to two dimensions, the choice becomes infinitely richer. Any analytic function of 
$ and X can serve as a superpotential, without spoiling the renormalizability of the two- 
dimensional model. We briefly discuss the issue of the continuous degeneracy of the soliton 
solutions in this setting (Sect. 2). 

The domain wall configurations, when viewed as depending on the Euclidean time r 
rather than on the spatial coordinate z, are known© to be the generating functions for 
amplitudes of multiple production of bosons at threshold by field operators. For the model 
discussed here these are the amplitudes for production of an arbitrary number of the (j> 
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bosons and k x of the x bosons at the corresponding thresholds, 

K k = (n^ x |0(O)|O) and A% k = <n^ x | X (0)|0) . 

In Sect. 4 we use the relation between the domain wall profile and the threshold production 
amplitudes in a twofold way: to point out a constraint on the solutions in the degenerate 
family stemming from the fact that they generate the same set of amplitudes and to find 
the multi-boson amplitudes at the tree level explicitly in the case of p = 4 where the explicit 
form of the solutions is available. 



2 Additional Integrals of Motion 

To begin with, we will discuss the occurrence of an additional integral of motion in the 
simplest example (|6]). In this model the creek equations have the form 

dW 

Let us introduce a "dual" function 

W($, X)=X- p (^-$ 2 —X 2 ) , (9) 



V A 2 r - , 

where 

A 

p = -. 

a 

The meaning of the word "dual" will become clear shortly. Equation @ assumes that the 
parameter p ^ 2. The case p = 2 is special and has to be treated separately. 
The dual function has the property 

dW dW 

wr^w, (10) 

where 0i j2 stands for {0, x}, £%j is the antisymmetric tensor of the second rank, and S is 
some scalar function. In the model at hand 

*(•.*) = ^ ■ (") 
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Now, if z is interpreted as "time", the dual function W is conserved along the trajectory. 
Indeed, 

~ dW ■ dWdW* 



For real solutions Eq. ( p!U|) implies that the right-hand side of Eq. (|T^) vanishes. 

Therefore, each particular real trajectory connecting the stationary points 1 and 2 of the 
superpotential is characterized by the value of the dual function on this trajectory. More 
exactly, dual functions conserved along the trajectory can be introduced for all solutions 
4>k(z) with the constant, i.e. z independent, phases of the fields <pk- By an appropriate 
redefinition of the fields we can obviously return to the real solutions. 

In the general case of non-minimal models the superpotential (restricted to real values 
of the superfields) is characterized by the gradient lines and the level lines. The latter 
correspond to fixed values of W. Two nets of lines - gradient and level - are locally orthogonal 
to each other. The level lines of the dual function are the gradient lines of the superpotential 
while the gradient lines of the dual function are the level lines of the superpotential. From 
this graphic interpretation it is intuitively clear that a dual function W must exist for every 
W, although, unlike the minimal two-field model, it is not always possible to find them 
analytically. The points where W (restricted to the real values of the scalar fields) develops 
maxima or minima are the singular points of W. The saddle points of W are the saddle 
points of W . 



2.1 Solitons in 1+1 dimensions 

If in four dimensions the choice of the superpotential is severely restricted by the require- 
ment of renormalizabilty (only polynomials which are at most cubic are allowed), in two 
dimensions any superpotential leads to a sensible quantum theory. If one takes a generalized 
Wess-Zumino model in four dimensions, with arbitrary number of fields and an arbitrary 
superpotential, and dimensionally reduces it to two dimensions, one arrives at N = 2 su- 
persymmetry in two dimensions. In two dimensions the domain walls become solitons - 
localized field configurations with finite energy. After quantization they are to be viewed 
as particles. The N = 2 supersymmetric theories with scalar superfields in two dimensions 
were extensively studied in Ref. 0. The case which is of most interest to us, a continuous 
family of degenerate solitons, seemingly was not discussed in this work. 

Since we are not limited now to polynomial superpotentials we can consider whole families 
of models. Consider for definiteness two-field models. Let us start from a certain model 
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with a superpotential W. One then can consider any other model with a superpotential 
Wnew = f(W) where / is an arbitrary function. If in the original model the dual function 
is known, it remains the same for the whole family. Indeed, Eq. (|10|) implies that 

— — = SijS-NEW — 7 W 5 'S'new = S (df/dW) 1 . (13) 

0<pi 0<pj 

In other words, W remains the integral of motion for real trajectories in any model belonging 
to the given family. 

The minimal two-field family, dimensionally reduced to D = 2, presents a simplest exam- 
ple where continuously degenerate soliton solutions exist. Now one can easily provide with a 
plethora of other interesting examples. For instance, one can consider superpotentials which 
are bounded from above and from below for real values of the superfields. In such models 
typically every soliton will appear as a member of a degenerate family of solitons, the degen- 
eracy being unrelated to any external symmetry. Generalizations of the sine-Gordon model 
fall into this category. Consider, for example, the superpotential 

W = -sin$ - sinX - a(sin$)(sinX) , (14) 

where a is a dimensional parameter, and all dimensional parameters are set equal to unity. 
This superpotential is periodic in $ and X; for a = it describes two decoupled fields (each 
of them presents a supergeneralization of the sine-Gordon model). If a ^ the fields <3> 
and X start interacting with each other. Inside the periodicity domain < $,X < 2tt the 
relief of the superpotential W is qualitatively similar to that of the minimal two-field model: 
—W has one maximum at$ = X = 7r/2, one minimum at $ = X = 3tt/2, and two saddle 
points at $ = tt/2, X = 3tt/2 and $ = 3tt/2, X = ir/2 (at least for small values of a). Any 
real trajectory (out of a continuous family of trajectories) starting at the maximum leads 
to a minimum. The only exceptions are two trajectories leading to the saddle points (the 
exceptional, or basic solitons). The masses of two exceptional solitons are 4 + 4a and 4 — 4a, 
respectively. Continuously degenerate solitons are bound states of two basic solitons, with 
mass 8, i.e. the binding energy exactly vanishes. 



Unlike the minimal two-field model all solutions in the model (14) have finite masses; 
there are no trajectories leading to abysses. 

The degeneracy is not lifted due to quantum corrections. This suggests that in every 
model with continuously degenerate solitons there should exist a dual description where 
the exceptional solitons (comprising the degenerate ones) appear as decoupled (i.e. not 
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interacting with each other) particles from the very beginning. This issue deserves further 
investigation. 



3 Solution for the Degenerate Walls 

Now we return to the minimal two- wall models (|B|). Our task is to find the family of the 
wall trajectories (i.e. <fi versus % for every given wall solution and every allowed value of W). 
Next we will find the wall solutions themselves, i.e. <p(z) and x( z )- 

3.1 General case: arbitrary ratio p = X/a 

It is convenient to introduce dimensionless field variables / and h as 

m , m . 

<P = ^rf, X = -7^ h , 
A v A a 

and to set the mass parameter m to one (it can be restored, if needed, from dimension). The 
BPS-saturation equations (f|) then take the form 

f = l-f-h 2 

dz 

— = — fh . (15 

az p 

By eliminating the variable z from these equations one finds the equation for the trajectory 
in the (/, h) plane, 



dh 2 fh 
The general solution of this equation can be written as 



(16) 



f = \-^--Ch\ (17) 
P * 

where C is an integration constant. It is connected with the integral of motion, Eq. (||]), by 
a simple proportionality relation, 

The full trajectory runs from the point (/, h) = (—1,0) in the "distant past" (i.e. z = — oo) 
to the point (/, h) = (1, 0) in the "distant future" (i.e. z = oo). The relation flT7f) determines 
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it piecewise: from (/, h) = (—1,0) to (/, h) = (0, /i max ) and from (0, h max ) to (1,0), where 
^max is the maximal amplitude of h on the trajectory (see Fig. 1). For the given value of C 
the value /t max is obviously determined from the solution of Eq. (|T7| ) for h with / being set 
to zero, 

T^max + C^max = 1 ■ (18) 



p-2 

Using the symmetry under h <-> — h it is sufficient to discuss only the trajectories with 
positive h. 

The freedom in C is in fact limited by the condition that the fields remain real along the 
whole trajectory. In this connection it should be noted that the BPS-saturation equations 
(0) are intrinsically non-analytic, thus the trajectory found for the real values of the fields 
cannot be continued to the the complex values of the fields. The requirement of the real 
trajectories translates into the requirement that h max is real positive and varies in the interval 
[0, 1]. An elementary inspection of Eq. fll8D yields the allowed domain for the constant C, 

2 

C* < C < +oo , C* = . 

2-p 

One can readily see that the trajectory with C = +00 has h(z) = and, thus, reduces 
to the well-known wall solution in the one-field theory, f(z) = tanh^. On the other hand 
the trajectory with C equal to the critical value C = C#, has h max = 1. Thus, it in fact 
describes two infinitely separated domain walls: one interpolating between (/, h) = (—1,0) 
and (0, 1) (the saddle point of W) and the other between (0, 1) and (1, 0). The configurations 
with the intermediate values of C interpolate between these two extremes; in a sense, they 
can be viewed as the solutions with the latter two walls at a finite separation (see Fig. 2). 
Remarkably, the degeneracy of the energy of the C > C* solutions implies that the latter 
two walls do not interact with each other. 

If p > 2 the point C = belongs to the allowed interval. For C = (i.e. W = 0) 
the trajectory takes an algebraic form. One can immediately find an explicit domain wall 
solution, see Eq. (0). As a matter of fact, it was obtained previously®. 

It should be noted that the value p = 2 presents a special case because of the singularity 
in Eq. (|17]) at p = 2. This singularity can easily be resolved, however, either by considering 
the limiting procedure in Eq. (|17|), or by solving the equation (|l^) separately for p = 2. The 
result is that the relation (|T7|) at p = 2 reduces to 

f = l + h 2 (\nh 2 -C) , (19) 
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1 

Figure 1: Few trajectories in the field space from the family of degenerate solutions (at p = 4). 
The numbers represent the values of the constant C for the corresponding trajectory. The 
heavy dots show the positions of the vacuum states. 



C=100 C=0 C=-0.9999 




Figure 2: The profile of the fields (solid) and x (dashed) as a function of the spatial 
coordinate z for three values of C (p = 4). It is seen that at C approaching the critical value 
= —1, the profile separates into two walls. 



10 



where C is another integration constant bound by the condition 1 < C < +00. 

The relation ( [T7| ) gives the trajectory of the field configuration in the (/, h) plane (see 
Fig. 1). In order to find the coordinate dependence of each of the fields, one can substitute 
/ as found from Eq. ( O ) into the second of the equations (|T^) and then obtain the solution 
in an implicit form, "in quadratures" 

P f dh 

z = -2lkM ' (20) 

3.2 Explicit solution for p = 4 

For arbitrary values of p there is no algebraic expression for the integral in Eq. ( pOf) in terms 
of known functions. The exceptional cases are 

1 2 

P = 2> 3> 6 > and 8 ' 
when the integral is expressed in terms of elliptic functions, and 

p = 1 and p — 4 , 

when the integral is elementary. 

The elliptic cases are rather cumbersome, while the case of p = 1 is trivial: at p = 1 the 
model considered here describes two fields <fi = (<f> + x)/ V% and \ — (x ~ 4>)/ V% that do not 
interact with each other. For these reasons we pursue here the explicit solution only for the 
exceptional case of p = 4. 

If we choose this specific value, p = 4, in Eqs. ([H]) and (p0|) , we readily find the explicit 
expressions for the fields / and h versus z, 

(u + iy + c (u + ly + c 

where u = exp(^ — zq) and zq is an arbitrary shift of the coordinate z. (Clearly, the solution 
can be centered at z = so that /(0) = 0, if z is chosen as 2 z = — ln(C + 1).) For p = 4 
the integration constant C is bounded by the condition 

-1 < C < +00. 

For completeness we also present the same solution in terms of the fields <p(z) and x{ z )i with 
the normalization factors restored (Fig. 2), 

e 2m(z-«,)_ C r_ 1 m 2 V2e m ^°)/ 2 . . 

A (e™(— o) + i) 2 + c AW A ^( e m( 2 - ao ) + 1 )2 + c 
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The constant C plays the role of a collective coordinate. The occurrence of another collective 
coordinate, z , is a trivial consequence of the fact that the wall solution spontaneously breaks 
the translational invariance of the original model in the z direction. 

At the same time, C is unrelated to spontaneous breaking of any symmetry of the model. 
It reflects the fact that the symmetry of the particular solutions at hand is higher than 
that of the model per se. As far as we know, Eq. (|22|) presents the first explicit example 
of a soliton family in the renormalizable 4 theory with a continuous degeneracy due to 
continuous deformations of the soliton structure. 

One can readily derive from the explicit solution above the limiting cases discussed in 
Sect. 2.1. At the critical value C — — 1 the solution in Eq. (pTT) degenerates into 



' = 1172 < ft2 = ^ (23) 

and describes the domain wall connecting the vacua (/, h) = (0,1) and (1,0). In order to 
get the wall connecting (/, h) = (—1, 0) with (0, 1) one can use the symmetry under / <-> — / 
and reverse the sign of / in the first equation in (^). In order to recover the one-field 
solution in the limit C +oo one has to accordingly adjust the coordinate shift zq as 
exp(— 2 zq) = C + 1. Then in the limit C — > +oo one gets h = and / = tanhz. 



4 Threshold multiparticle amplitudes 

In this section we will take advantage of the explicit wall family solution found above in order 
to extract certain predictions for the high-order behavior of the multiparticle amplitudes at 
thresholds. The corresponding analysis for the one-field Wess-Zumino model was carried out 
in Ref. J|]. In the one-field model there is little distinction with the non-supersymetric case 



(for a review see Ref. [0). The two-field model is much more interesting since it reveals a 
new pattern. 

In what follows we will need to know that in the vacuum with = m/\ and x — (from 
which our wall trajectory starts) the mass of the quantum is equal to 2m while the mass 
of the x quantum is equal to 2m/ p. In the case p = 4 to be analyzed below the mass of the 
X quantum is equal to m/2. The same is valid for the vacuum = m/A, x = 0, from which 
the trajectories originate. 
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4.1 An overview of the formalism 



The solutions to the field equations, in particular the domain wall solutions, are directly 
related to multiparticle amplitudes, by virtue of the formalism developed by BrownGH (for 
a more recent review see Ref . ||14|| ) . Being adapted to the present problem of two fields the 
formalism is constructed as follows. Consider for definiteness the amplitude (n, fc|0(O)|O) 
describing the production by the field operator <f)(x) of a multiparticle state consisting of n 
on-shell bosons of the field with 4- momenta p a (a = 1, . . . ,n) and k on-shell bosons of the 
field x with 4-momenta p b (b = 1, . . . , k) in a vacuum |0) of the theory. According to the 
standard reduction formula this amplitude is expressed through the response of the system 
to external sources p^x) and p x (x), coupled to the corresponding fields as p ( / > (fi + p x x i n the 
Lagrangian, 



(n, A#(0)|0> 



n . 

n / 

2=1 Pl^ m i J 



n 2 iim 9 / ^ e 

b=1 Pt^ml J 



ipbXb (ml 



(fx a e ipaXa (m 2 6 
-Pb) 



Pi) 



X 



5p x (x b ) 



Sp4,(x a ) 

(Oout\m\Oin) P *' P % = o, Px = 



(24) 



where and m x are the masses of the respective bosons in the vacuum considered. 

The classical response, i.e. the classical solution of the field equations in the presence of 
the sources generates, through Eq. (|24j), the tree-level amplitudes, which we will be mainly 
concerned with here. Moreover, as will be seen, the configurations, of the type of the domain 
walls, depending on only one variable, are related to multiparticle production exactly at the 
threshold, i.e. at the spatial momenta of the produced particles exactly equal to zero. In 
this situation it is sufficient to consider the response of the system in Eq. (^4j) to spatially 
uniform time-dependent sources, 



Prj) {t) = p^) e*** , p x (t) = Px (uj x ) 



and take in the very end the on-shell limit in Eq. 



by tending oj^ to and uj^ 



The spatial integrals in Eq. (|24|) then give the normalization spatial volume, conventionally 
set to one, while the time dependence with the fixed functional form of the sources implies 
that the propagator factors and the functional derivatives enter in the combination 



Pi) 



5 



I 2 2- 



5pj,(x a ) 
6_ 

Sp x {x b ) 



(m. 



S 



5 



5 



5a<p(t) 
5 

5a x (t) 



(25) 
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where 

am = 4^, ax{t) = M^± (26) 

mi — ie — ui — it — uj^ 

coincide with the response of free fields to the external sources. For finite amplitudes of 
the sources the response is singular in the on-shell limit — *> m^, u x — > m x . Therefore, 
following BrownO, the amplitudes of the sources should be taken to zero in this limit, so 
that the factors a^t) and a x (t) are finite, 

tty(t) -> a e im ^ , a x (t) -> a x e im ^ . 

Thus, for the purpose of calculating the multiparticle amplitudes at the tree level one 
looks for a solution of the classical field equations with no sources. The only information 
about the sources left from the above limiting procedure is that the sources drive only positive 
frequencies in the fields, thus the condition for the sought solution is that it should contain 
only the positive frequency part with all harmonics being multiples of e* m< *' and e tnixt . The 
latter condition is equivalent to requiring that the solution for the fields goes to the classical 
vacuum at infinity in the Euclidean time r = Imt — > +00. The multiparticle amplitudes are 
then given by the derivatives of the solution <f)(t), 

<- s <"■ * wo)|0> = fe^))" (aSij) * m ° ' (27) 

Since the equations for the fields and x are coupled, one simultaneously finds the solution 
for the field x an d, thus, the amplitudes for the multi-boson production by x, 

A ik s <„, %«»io) = " (^y 1 m \ at __^__ . m 

The operational procedure for calculating the amplitudes is therefore as follows. First, 
one obtains the solution of the Euclidean classical field equations depending only on time r 
and approaching, at r — > +00, the vacuum state (0o, Xo) m which the scattering theory is 
considered. Then, the solution is expanded in the harmonics of e~ m <> > r and e" 



E F n , k e-( nm * +km ^ T , X (r)= E H n , k e~( nm * +km ^ T , (29) 



n=0, k=0 n=0, fc=0 

where F n ^ and H n k are the coefficients of the expansion. 
Note that 

^0,0 = 0o and H = Xo , 
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while the coefficients of the appropriate first harmonics are identified as the described above 
factors a«£ and a x , 

Fi, o = d(p , Hq > i = a x . 
Then, according to Eqs. ( p7f) and (p8|) , the amplitudes are expressed as 

<k = n\kl-J^- , Ax h = n\kl^£- . (30) 



Before closing this discussion of the general formalism we would like to emphasize that 
the latter equations can be also viewed as a constraint on any solution approaching a vacuum 
state (00) Xo) at r ^ +oo. Namely, if up to the appropriate linear terms the fields behave 

as 

0( r ) = o + ae- m ^ + ... , X (r) = X o + be™** + . . . , 

then the subsequent terms in the expansion of the fields are fully determined in terms of a, 
b and the fixed set of the amplitudes and A x , 

At k a n b k Al k a n b k , 

jp _ n , K TT n,K (q-i\ 

r n ,k — rrj , ti nk — — . [61) 

Furthermore, the absolute normalization of the coefficients a and b is rather a matter of 
convention. Indeed, under a shift of r, 

t -> r -t 

these coefficients change as 



ae m * T \ b^be m * 



TO 



Thus, the only parameter that distinguishes between essentially different solutions approach- 
ing the same vacuum state is the ratio 



Therefore, in a general two-field theory a family of solutions approaching a vacuum state 
at r — > +oo is parametrized by a single parameter c. This parameter is in one-to-one 
correspondence with the integration constant C, or the value of the dual function W on 
the trajectory. The coefficients of the expansion (EH) are then fixed by the multiparticle 
amplitudes. 
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4.2 Multiparticle amplitudes in the supersymmetric model 

The domain wall solutions discussed in Sect. 2 can be directly applied to calculating the mul- 
tiparticle amplitudes. To this end one should consider the fields depending on the Euclidean 
time t rather than on the spatial coordinate z. Since this amounts to a trivial relabelling 
of the variable, we retain here the notation z for the variable. We also use the notation 
a = Fi q and b = Hq \. Every BPS-saturated solution from the family under consideration 
approaches at z — > +00 the vacuum at (0, x) = ( m /% 0). Thus, this is the vacuum state 
in which the multiparticle amplitudes are generated by the solutions. Remember that the 
masses of the particles in this vacuum are expressed in terms of the parameters of the model 
as = 2 m, m x = 2 (a/X)m = m^/p. 

In the case of arbitrary ratio p one can use Eq. (|T7| ) to obtain a relation between the 
coefficients a, b and the constant C . Indeed, at z — > +00 the field x( z ) g° es to zero as 

X (z)=be™** + ..., 
corresponding in the dimensionless variables to 

h(z) = y^e-^. 
m 

The linear in e~ 2z harmonics in f(z) arises from the term in Eq. ( |T?D with the constant 
C, from where one finds the coefficient a of the linear in e~ m 4> z harmonics in (f>(z) ((f>(z) = 
f + ae~ m * z + ...) as 

Cm/ \f~\a b\ P 



, ■ (32) 

In connection with the derivation of the latter relation it should be noted that for p > 2 
the harmonics with exp(— m^z) is not the leading one in the field <f>{z) at large z because 
of the presence of the second harmonics with the mass of x~- -^0,2 exp(— 2m x z). Thus, the 
presented derivation, strictly speaking, is justified only at p < 2. However, the relation ( |52"D 
is also applicable at p > 2 since the fields and the coefficients of their expansion are analytic 
functions of the couplings, and the relation ( |32"1) can be analytically continued from p < 2 to 
the domain p > 2. 

It can be also noted that the solution in Eq. (|7]) with C = has a = 0, according to 
Eq. (j32|). Thus, it generates only the amplitudes of multiple production of the bosons of 
the field x by either the operator x(0) or (f)(0). For this reason it expands in the harmonics 
determined only by the mass of x- 
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Furthermore, for a rational ratio p the masses m$ and m x are also in a rational proportion. 
Thus if only the z dependence of the fields were known, there would be an ambiguity, at least 
in some harmonics, in separation between the production of the x bosons and the bosons. 
However, this ambiguity is resolved if the dependence of the solution on C is known, by 
using Eq. Q52|), which shows that the constant C serves as a "tag" for a boson. The power 
of C in the given harmonics gives the number of the </> bosons in the amplitude generated 
by this harmonics. 

We illustrate this method for our explicit solution in the case of p = 4 and we also find 
explicitly the amplitudes k and A\ k in this case. 

Setting for definiteness Zq = in the explicit solution in Eq. (|22|), we find the coefficient 
b determining the rate of approach of the field xi z ) to its vacuum value (zero), 

b = — - — . (33) 
Furthermore, using the relation (|32| ) at p = 4 we also obtain the coefficient a, 

77? 

a = -2C- . (34) 
A 



We then expand the expressions for the fields in Eq. ( |22|) in powers of C and, finally, each 
term of this expansion in powers of e~ mz . In this way we get 



a „-(2n.+ |) mz 



<K*)= E Fn, k C n e 

n=0, k=0 
oo 

X(z)= £ H n , k C n e-( 2n+ -*) mz , (35) 

n=0,fe=0 

where the coefficients T n ^ (H n ,k) are non-zero for even (odd) k, as is expected from the 
symmetry of the model at hand. The explicit expressions for these coefficients are 

X { 1} (2n)!(*/2)! ' ( + } ' 
_ 2V2m +(& _ 1)/2 r(n + l/2) [2n + (fc - l)/2]! 
^ n ' fc " A 1 j v^Fn!(2n)! [(Jfc - 1)/2]! ' 1 j 

The combinations T n , u C n and 7i nj k C n are identified as respectively the coefficients F n> k 
and H n k in the general expansion of Eq. (p9|) . The latter coefficients are related to the 
multiparticle amplitudes as given by Eq. fl3lD . Using the explicit expressions in Eqs. ( p3|) 
and (^4|) for 6 and a, one arrives at the relation between the amplitudes and the found 



17 



x(pi) 

x(Ps) 

x{pd 



4> + (k-i) x 



Figure 3: A graph with four \ bosons originating from a single virtual <fr, which is singular at 
threshold if m ( f ) /m x = 4 unless the scattering amplitude a (filled circle) vanishes at threshold. 
Both the open and the filled circles represent the sum of tree graphs. 




coefficients JF n fc and H n ,k- Namely, 

n+k / \ \ n+k—1 



n, k 



-lYn\k\ f AY = f-D fc / 2 [AY n\k\(2n + k/2)l 



Al u = (-l) n n\k\ 



mj 2"+ 3fc / 2 v ' \mj 2 n+3k / 2 - 1 {2n)\{k/2)\ 
\ rt n k 



n,k v ±j y m j 2n+3k/2 

= r _ nM /2 /AV^ 1 k]T (n + l/2) [2n + (fc-l)/2]! 
1 ' \m) v ^F2«+ 3 ( fc - 1 )/ 2 (2n)! [(k - l)/2]! ' 1 ' 

This concludes our calculation of the threshold amplitudes in the minimal two- filed model. 



4.3 New zeros 

A remarkable property of our result is that the amplitudes are finite even though a state of 
four x bosons at threshold is degenerate in energy with one <fi boson. In other words, any 
graph, where four final % bosons with the four-momenta pi, p 2 , P3 and p<± originate from a 
single line of (see Fig. 3), contains the factor 

oi, 4 (pi, P2, Pa, Pa) , qq x 

7 To 9 • l^°J 

(Pi +P2 + P3 +PaY ~ m% 

Here ai^fyi, p 2 , P3, p±) is the conventional Feynman scattering amplitude for the process 
<f) — > 4x- At threshold the denominator in Eq. (|38|) goes to zero, and the graph becomes 
singular unless the Feynman amplitude ai j4 (pi, p 2 , P3, p±) also vanishes when all four mo- 
menta are at threshold. The latter cancellation indeed takes place in the model considered 
here. This can be seen by examining the amplitude Af 4 , which exactly is the threshold limit 
of the expression in Eq. (j38|). Indeed, the amplitudes A^ and A x considered here are the 
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matrix elements of the field operators in Eqs. fl27|) and (p8|). These matrix elements have 
the propagators of the final on-shell bosons amputated, but the propagator of the incoming 
virtual field is not amputated and remains included in the corresponding amplitude A. The 
conventional Feynman scattering amplitude is thus obtained by multiplying the amplitude 
A by the inverse propagator of the incoming line. In the case of the process — > 4% the 
inverse propagator of the incoming <fi is vanishing at threshold of four \. Thus, the Feyn- 
man scattering amplitude also vanishes. Clearly, this cancellation can be also verified by an 
explicit calculation of the tree Feynman graphs. 

This cancellation can be extended to a general case of an arbitrary even integer value 
of p, with the exception of p = 2. In the latter case the exponential behavior of % at 
z — > +oo generates a non-exponential dependence of through the logarithm in Eq. (fl9|), 
which implies that in this case a resonance between the degenerate states does take place. 
For all other values of p the coefficients in the expansion of the type as in Eq. ( p9"D can be 
constructed by iterations and are non-singular. This means that for the values of p where 
the resonance could potentially occur, i.e. — > p x for even integer p, it actually does not 
take place due to vanishing of the corresponding Feynman scattering amplitude. 

5 Conclusions 

In summary, in a class of supersymmetric models with the continuously degenerate family 
of BPS domain walls (with real trajectories) an additional integral of motion is observed. 
The occurrence of this integral allowed us to find a generic solution from the family in 
quadratures, while for a specific ratio of the coupling constants the whole wall family is 
obtained in the closed form in terms of elementary functions. We then further utilize the 
result for deriving the multiparticle amplitudes at threshold in the minimal two-field Wess- 
Zumino model. The threshold amplitudes are calculated in a closed form for p = 4. In the 
course of the calculation we have found an unexpected cancellation of the tree graphs for the 
Feynman amplitude of the process <fi — > 4 \ at the threshold, due to which cancellation the 
multiparticle threshold amplitudes are finite. We also conclude that the same cancellation 
takes place for the process <fi — > px at arbitrary even integer p, except p = 2. The relation of 
this cancellation to additional integrals of motion is yet to be studied. It can be also noted 
that the nullification of the amplitudes is somewhat reminiscent of the general property of 
nullification!!!!' for the on-shell processes 2 — > many at the threshold in scalar theories. 
Interesting phenomena occur when the models at hand are dimensionally reduced to 
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D=2. The two-dimensional theories thus obtained have extended supersymmetry, N = 
2. A continuous degeneracy of the soliton family (persisting with all quantum corrections 
included) reflecting the possibility of the continuous deformations of the solution profile can 
be seemingly interpreted in this case as the existence of decoupled "basic solitons" . Revealing 
these decoupled basic solitons in an explicit form and studying their properties is an obvious 
next problem to be dealt with in the given range of questions. 
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